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We present several strategies for searching for supersymmetry in dijet channels that do not ex- 
plicitly invoke missing energy. Preliminary investigations suggest that signal-to-background ratios 
of at least 4-5 should be achievable at the LHC, with discovery possible for squarks as heavy as ~ 
1.7 TeV. 
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I. INTRODUCTION 

The LHC is set to explore the physics of the weak 
scale, whatever it should turn out to be. Supersymme- 
try is one of the leading candidates and enormous effort 
has been dedicated to studying missing energy signals 
that characterize almost any weak-scale supersymmetric 
model. However, supersymmetry searches will be chal- 
lenging and disentangling the supersymmetry parameters 
will be more difficult still. 

In light of the above, it is imperative to study every 
possible channel in order to optimize our chances of dis- 
covering new physics and understanding the underlying 
theory. In this regard, events with the lowest multiplicity 
may be the simplest ones with which to make headway 
on the inverse problem. 

Although two-jet events with missing energy have been 
studied at the Tevatron [l[, they have been less promi- 
nent in LHC studies. ATLAS has shown that two jet 
events can be useful for certain SUSY models, both for 
discovery and for constraining superpartncr masses 0, 
but recent ATLAS and CMS studies have focused more 
heavily on the more challenging cascade decays. In this 
paper we study one novel and two existing kinematic vari- 
ables that can be used to capture dijet missing-energy 
events without explicit reference to missing transverse 
energy. We find that pairs of these variables can be used 
to give signal-to-background of at least 4-5, indicating 
that these variables are worth exploring with a full de- 
tector simulation^. 

Dijet events are worthy of attention as a potentially 
clear window into parameter space. They are not results 
of complicated cascade decays but arise simply from two 
squarks decaying to two quarks and two neutralinos. Be- 
cause we know the identity of the particles involved and 
because there are so few, the signal is relatively straight- 
forward to interpret. For example, with sufficient inte- 
grated luminosity, these events alone can be used to con- 
strain the squark and neutralino masses. Dijet studies 
along the lines explored here may usefully supplement 



1 Such a study has been started by CMS Q. In addition, ATLAS 
is currently engaged in an updated dijet study 



recent analyses dedicated to distinguishing SUSY from 
other models using events with at least three jets Q. 

The kinematic variables we consider are constructed 
from the two jets' momenta alone. These variables 
should have different systematic uncertainties than miss- 
ing transverse energy since they pick out slightly different 
events and are based on different measurements. At the 
very least, then, the searches we suggest should be worth- 
while as cross-checks of standard searches. The variables 
we use may also be useful for optimization when signal- 
to-background is relatively low. 

The searches we describe will be most effective when 
squarks are pair-produced in abundance and have large 
branching ratios to decay directly to the lightest neu- 
tralino, which requires that squarks are lighter than the 
gluino so that cascade decays through gluinos arc ab- 
sent. Because t-channel gluino exchange is an important 
source of squark pair production, the lighter the gluino 
the more prominent the signal. For the parameter points 
considered below, we find the signal is cut by a factor of 
~ 6 — 7 when the gluino decouples. Fortunately, compa- 
rable gluino and squark masses are a feature of a large 
class of models - most notably high-scale models where 
the heavier gluino mass feeds into the squark mass. We 
focus on such models in this study. 



II. ANALYSIS DETAILS 

Before getting to the dijet properties that will be the 
focus of our study, we consider the effectiveness of I^t 
and , the missing transverse energy obtained from 
the dijet system alone. After requiring the sum of the 
two jets' pt's to be greater than 500 GeV, event rates 
and signal-to-background ratios for one particular SUSY 
point are presented in Table U (details regarding event 
generation and cuts arc given below). Neither variable 
suffices for a clean search, but we observe that the S/B 
values obtained using p^T3i.re essentially identical to those 
obtained using I^t ■ This analysis suggests that, in the 
two-jet channel at high pT, nothing is to be gained by 
using full rather than kinematic variables associated 
with the two jets alone. 

We now present three dijet variables that can be used 
to separate signal and bacground, with ~ 1% of signal 
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TABLE I: For dijet events passing the cuts described in the 
text, the dependence of the signal cross section and signal- 
to-background (SjS) on a variable IpT cut (top), and on a 
variable _^TCut (bottom). All energies are in GeV. 

events passing all cuts. 

• a: which we define as the ratio of the of the 
second hardest jet and the invariant mass formed 
from the two hardest jets, 



As far as we know, this variable has not been con- 
sidered previously. Background events generally 
trail off at a = 0.5, whereas supersymmetry events 
with invisible decay products can easily have larger 
a. Large a tends to arise in events in which the jets 
are not back-to-back. As one extreme example, if 
the two jets are nearly aligned, their invariant mass 
can be quite small, leading to very large a. 

Because of the background's sharp drop-off around 
a = 0.5, this variable is potentially useful as a diag- 
nostic tool for analyzing two jet events and cleanly 
separating signal events from QCD. 

• A(/): the azimuthal angle between the two hardest 
jets. Azimuthal angle is often used in conjunction 
with missing transverse energy, and A(/) was among 
the variables used in the diiet SUSY search at DO 

Q. 

• Mt2 H: which is defined for events in which two 
particles of the same mass undergo identical semi- 
invisible decays, as 

Mt2{x) = , min {max[mT(pi, ^i, x), mT(p2, ^2, x)]}, 
ni+n2=PT 

(2) 

where pi and p2 are the momenta of the visible par- 
ticles, is the missing transverse momentum of 
the event, and tot is the transverse mass function, 
which depends on an assumed value x of the invis- 
ible particle's mass. In calculating Mt2(x) use 
the missing transverse momentum as determined 
by the dijet system alone. 

If X is taken to be equal to the mass of the in- 
visible particle, the Mt2 distribution will have an 
endpoint at the mass of the decaying particle. Not 
knowing this mass, Mt2 endpoints still constrain 
the masses of the decaying and invisible particles, 
as emphasized in [l^l and used below. 



We consider these variables singly and in tandem. 
We find the first two variables are useful in that one 
can choose parameter-independent cuts that give sizable 
S/B, whereas the last variable, though more parameter- 
dependent in its optimization, might ultimately maxi- 
mize S/B. Since the advantage is not overwhelming, we 
expect all the variables could prove useful, either at the 
trigger or analysis level. Because they are dimcnsionless, 
the first two variables might have the further advantage 
of being less sensitive to absolute energy scale, and might 
therefore have lower systematic errors. 

For all our analyses, we select events in which exactly 
two jets have pr > 50 GeV, with no isolated leptons, 
photons, or r jets. One could attempt to achieve better 
background rejection by an additional veto on extra jets 
with lower px- In general, we have chosen felicitous cuts 
but have not pursued a careful optimization, which will 
be more appropriate at the full-detector-simulation level. 

A gluino that is only slightly heavier than the squarks 
arises naturally in models with supersymmetry broken at 
a high scale, as renomalization-group effects prevent the 
squarks from being hierarchically lighter than the gluino. 
For our analyses we specify parameters at the high scale 
and use the SUSY-HIT package 0] to calculate super- 
partner masses and decay branching ratios. In the rel- 
evant parameter regions, the signal depends strongly on 
Ml/2, the unified gaugino mass at the high scale, and is 
less sensitive to Mq, the unified scalar mass, because the 
squark mass is dominated by gauge-loop contributions. 
We set the other SUSY parameters to be tan/3 = 10, 
Aq = 0, and /i > 0. 

The backgrounds included in our analyses are QCD, 
{W — !■ lv)/{Z z^I7)-|-jets, and tt. We have checked 
that diboson-f jets production does not significantly mod- 
ify our results. The QCD and tt samples were generated 
with Pythia 6.4 0], and Z/WViets with Alpgen 2.12 
0. Fully showered and hadronized events were then 
passed to the PCS 4.0 detector simulator 0, with the 
energy smearing in the hadronic calorimeter given by 
IS.E/E ~ O.S/i/S/GeV and the calorimeter granular- 
ity set to (A0 X A77) = (0.1 X 0.1). Jets were defined 
using a cone algorithm with Ai? = 0.4. 

A K-iacto^ of 2 is applied to the QCD sample, but no 
isT-factor is used for W/Z production, because the most 
important contributions come from W/ Z-\-2 jets, which 
arc not enhanced at NLO [13. (After cuts, W/Z produc- 
tion ends up being the dominant background to SUSY 
dijet events, so to include a A'-factor one can simply di- 
vide our signal-to-background ratios by K.) For tt we 
use cr = 830 pb as the NLO production cross section [ll| ■ 
Including the K factors our samples sizes are ^ 0.8 fb~^ 
for QCD, ~ 20 fb"! for tt, and ~ 100 fb'^ for W/Z. Ap- 
propriate generator-level kinematic cuts were imposed to 
obtain the QCD and W/Z samples. 

SUSY samples were also generated with Pythia. For 
each parameter point we use Prospino 2.0 [l3| to calculate 
an appropriate i^-factor from the NLO cross section for 
squark pair production [l3| . 
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FIG. 1: SUSY, Z — » !^!/+jets and QCD rates for passing the cuts described in the text, as functions of a (left), A0 (middle), 
and AfT2(0) (right). 



III. RESULTS 

The plots in Figure [T] suggest that appropriate cuts 
on a, A0, and/or Mti can suppress both the QCD 
background and the dominant background after cuts, 
(Z — *■ w)+jets. The SUSY parameter point used here is 
(Ml/2, ^o) = (300, 100) GeV, and we impose a hard cut 
on the sum of the two hard jets' transverse momenta, 

VTx + PT2 > 500 GeV. (3) 

To streamline the analysis, events were required to have 
JfjT > 100 GeV for Figure, [Hand at least one of a > 0.5, 
A(j) < 27r/3, and > 100 GeV for Figured Remov- 
ing these requirements does not affect the results once 
optimal cuts on a, Acj), and/or Mt2 are made. 

Evidently signal dominates over background for 
a > 0.5, A(/)< 27r/3, and Mt2 > 300 GeV. We will soon 
see that a, A0, and Mt2 can be used to discriminate 
signal from background by themselves, but first we point 
out that cuts on these variables can improve an anal- 
ysis based on or ^It ■ For example, the combi- 
nation (a > 0.45, fir > 300 GeV) selects 315 signal 
events per fb^^, with S/B = 4.3. The combination 
(A0 < 27r/3, > 450 GeV) gives a somewhat lower 
S/B (3.1), but with more events (429). An Mt2 cut 
of 450 GeV gives the largest S/B of all (5.0, with 304 
events), and in fact there appears to be no benefit in 
supplementing the Mt2 cut with the ^It cut. 

Figure [2] suggests that each of a, A0, and Mt2 can 
be used independently to observe a clear signal, without 
employing ^^at all. Well-chosen cuts give ^ a few x 10^ 
signal events after 1 fb~^, with S/B 3 — 5. 

Figure[2]also shows how the three variables can be used 
in pairs to improve S/B in conjunction with the signal 
event-rate. We again find that Mt2 seems to dominate 
a little, but since we do not know if this is the cleanest 
variable to use in practice, which can be determined nly 
after a full detector simulation, we present all combina- 
tions. Any two on their own can potentially give a robust 
signal. 

As an exanple, we consider the combination A0 < 
27r/3 and a < 0.45, which gives a good S/B and a decent 
event rate. As stated earlier, we do not optimize cuts, 
but we use this combination that works rather well. 
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Mt2{0) (GcV) 

FIG. 2: For events preselected as described in the text, the 
dependence of the signal cross section and S/B on a variable 
a cut (top), a variable cut (middle), and a variable Mt2 
cut (bottom). 



With those cuts in place. Figure [3] shows signal and 
background events binned in the sum of the two hard- 
est jets' transverse momenta. We see that Z-|-jets is the 
dominant background, followed by VF-l-jets. A total of 
four QCD events with pij, + P2t < 500 GeV passed the 
cuts, out of a sample corresponding to over 1.5 fb~^ of 
integrated luminosity, divided by the K factor. A higher 
luminosity sample would be needed to get a better esti- 
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FIG. 3: Signal and background rates after the cuts A(j> < 2n/3 
and a > 0.45. The QCD background is not included for 
P1T+P2T < 500 GeV. We take (Af^a, Afo) = (300, 100) GeV. 



FIG. 4: The mT2 distribution for signal and background, af- 
ter the cuts described in the text. We take {M1/2, Mq) = 
(300, 100) GeV. 
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TABLE II: The efficiencies e for signal events to pass the 
successive cuts, taking (Mi/2,Mo) = (300, 100) GeV. 



mate of the QCD background, but it seems safe to say 
that the W and Z backgrounds are more important. 

In Figure [3] we see that S/B is cleanest at high pr- 
Of course the optimal px cut depends on underlying pa- 
rameters that are not known a priori, but a scan at high 
Pt should help maximize S/B. For the chosen param- 
eter point, cutting above + = 550 GeV gives 
S/B = 4.9, with an average of 205 signal events after 1 
fb"""^ . Table lU shows the efficiencies with which the SUSY 
events pass the successive jet multiplicity, pxi +PT21 
and a cuts. 

The final efficiency is lower than that for SUSY 
searches with additional jets, and so despite the differ- 
ent systematics SUSY might well be discovered in other 
channels first. Moreover, the dijet channel is relevant 
only for certain models. On the other hand, this analysis 
picks out particularly simple events- two squarks decay 
to produce two jets and two neutralinos. If these events 
do occur it would certainly be worthwhile to study them 
in isolation. 

For example, with enough luminosity these events 
alone can be used to obtain a simple constraint on the 
squark and neutralino masses, using the Mt2 event func- 
tion [3] introduced above. If one can ignore all visible 
particles in the event except those in the two jets, one 
expects the endpoint 



MT2{0)r. 



771? - 77l|o 



(4) 



For the parameter point under study, the predicted end- 



point turns out to be 619 GeV if we use the mass of 
the right-handed squarks, which are the ones that decay 
predominantly to XiQ- Figure 2] shows the 777.^2(0) dis- 
tribution for 10/fb of data, with the cuts of Table HIl 
imposed. A sharp drop-off leading up to ~ 620 GeV is 
evident, consistent with expectations. The spill-over to 
larger values is mostly due to the effects of extra jets not 
included in the calculation of the missing transverse en- 
ergy (in calculating Mt2 (0) we use the missing transverse 
momentum as determined by the dijet system alone). 

The (a, A(/)) analysis we have described can be effective 
for higher-mass searches as well, with the cut on the sum 
of the two jets' transverse momenta increased appropri- 
ately. Table Hill gives results for other parameter points, 
with the cuts on pij^ + P2t again chosen to give robust 
values of S/B. The Mq values are chosen to be near the 
lower bounds below which a stau LSP results. Provided 
that the squarks remain lighter than the gluino, increas- 
ing Mq lowers the event rate somewhat but not dramat- 
ically. For {Ml/2, Mo) = (300,300) GeV, for example, 
the same cuts used for the (M1/2, Mq) = (300, 100) GeV 
point give 195 events after 1 fb"^, with S/B = 4.7. 

Taking S/^/B > 5 as the relevant criterion, our results 
suggest that discovery through the dijet channel should 
be possible for squark masses up t abut ~ 1700 GeV after 
100 fb~^ of integrated luminosity. By the same measure, 
discovery for lighter squark masses, ^ 600 GeV, should 
be possible after ~ a few x 10^ pb~^ or less. It may be op- 
timistic to focus on S/ ^/B as a discovery criterion, as do- 
ing so assumes that the background is fully understood. 
However, it is worth pointing out that (1) events with 
leptonic Z decays will provide some experimental handle 
on the dominant background, Z-|-jets, and (2) the shapes 
of the PT1+PT2 distributions for signal and background 
events passing the a and A(p cuts are quite different (see 
Figure . The excesses obtained in our analysis would 
lead to a prominent bump in the measured distribution, 
which would not be accommodated simply by rescaling 



5 



(Ml/2, Ma) 


("lg,"l5„) 


J]pt cut 






^sy (fb) 


S/B 


(300, 100) 


(716, 640) 


550 


9.5 X 


10" 


3 




205. 


4.9 


(450, 100) 


(1040, 918) 


800 


7.9 X 


10" 


3 




21.3 


4.7 


(600, 150) 


(1358, 1195) 


1050 


8.1 X 


10" 


3 




4.07 


5.0 


(750, 200) 


(1669, 1465) 


1250 


9.6 X 


10" 


3 




1.17 


4.8 


(900, 200) 


(1965, 1726) 


1450 


1.0 X 


10" 


2 




0.37 


3.5 



TABLE III: Efficiencies, event rates, and signai-to- 
bacfcground ratios for various SUSY parameters, using the 
cuts described in tlie text. Ail masses are in GeV. 

the background. 

IV. CONCLUSIONS 

We have studied several kinematic variables that can 
be used for dijet SUSY searches, and found that they 



give reasonable signal-to-background ratios. Dijet events 
can be used to constrain SUSY mass parameters should 
the type of supersymmetry model we have considered be 
correct. Studies of Z+jet events with leptonic Z decays 
will give a better understanding of the background and a 
more reliable extraction of signal from background. For 
the future, it would be useful to see how well the lessons 
here can be applied to develop multijet searches that do 
not rely on full missing energy. 
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